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ABSTRACT 



We present a new model for the full shape of large-scale the power spectrum 
based on renormalized perturbation theory. To test the validity of this prescription, 
we compare this model against power spectra measured in a suite of 50 large volume, 
moderate resolution N-body simulations. Our results indicate that this simple model 
provides an accurate description of the full shape of the power spectrum taking into 
account the effects of non-linear evolution, redshift-space distortions and halo bias for 
scales k < 0.15 /iMpc^^, making it a valuable tool for the analysis of forthcoming 
galaxy surveys. Even though its application is restricted to large scales, this prescrip- 
tion can provide tighter constraints on the dark energy equation of state parameter 
wbe than those obtained by modelling the baryonic acoustic oscillations signal only, 
where the information of the broad-band shape of the power spectrum is discarded. 
Our model is able to provide constraints comparable to those obtained by applying 
a similar model to the full shape of the correlation function, which is affected by dif- 
ferent systematics. Hence, with accurate modelling of the power spectrum, the same 
cosmological information can be extracted from both statistics. 

Key Vifords: large-scale structure of Universe - cosmology: theory. 



1 INTRODUCTION 

Over the past decade, new cosmological observations have 
revolutionised our view of the Universe by showing increas- 
ing evidence that it is currently undergoing a phase of ac- 
celerated expansi on. The Hubble diagram of Type la su- 
pernovae (SNlal (iRiess et all 1 19981 : iPerlmutter et al.l 1 19991 : 
iKowalski et al.ll2008l ). together with other independent data 
sets like observations of the temperature fluctu ations in 
the cosmic microwave background (CMB) (Hinshaw et al.l 
l2003l : ISpergel et ai]|200l l2007l : iKomatsu et al.ll2009D conT 
bined with measurements of the large scale structure 
(LSS) of the galaxy distribution (jEfstathiou et al. 20021: 
Percival et al. 2002; Tegmark et al. 2004; 'S anchez et al.l 
2006 . 2009 : Percival et al. 2009 : Reid et al. 2oTdi ]. have es- 



tablished the picture of a geometrically flat Universe whose 
energy content is composed of about 4% baryonic mat- 
ter, 21% of cold dark matter and 75% of a mysterious 
component called Dark Energy (DE), responsible for the 
current acceleration of the expansion rate. The nature of 
dark energy is one of the most important open questions 
in modern cosmology as the comprehension of its nature 
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has deep implications for our understanding of fundamen- 
tal physics. In the past few years a substantial efi'ort has 
been devoted to constrain the properties of dark energy 
as accurately as possible in order to distinguish between 
possible alternative descriptions o f this phenomenon (see 
ICopeland. Sami fc Tsuiikawal l200q . for a review). This ef- 
fort will continue in the coming years with the construc- 
tion of the next generation of large galaxy surveys, like 
the Panoramic Surv ey Telescope fc R apid Response Sys- 
tem (Pan-ST ARRS, iKai ser et al.l 120021 ') . the Dark Energy 
Survey (DES. lAbbott et al. 2005 ), the Baryonic Oscillation 
Spectroscopic Survey (BOSS. ISchlegel. White fc EisensteinI 
[20091), the Hobby Eberly Te lescope Dark Energy Experiment 
(HETDEX. iHill et a]ir2004l ') and, on a longer tim e scale, the 
space based Euclid mission (|Cimatti et al.ll2008r ) . These sur- 
veys will cover volumes much larger than current datasets, 
providing a dramatic improvement in the accuracy of the 
constraints on the values of cosmological parameters, in par- 
ticular the dark energy equation of state parameter deflned 
as u;de = Pde/pde, where pde and puE denote the pressure 
and the energy density of this component. 

These new surveys will provide a clear picture of 
the Baryonic Acoustic Oscillations (BAO), a signature im- 
printed in the large scale galaxy distribution by the acous- 
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tic fluctuations in the baryon-photon fluid prior to recom- 
bination. These can be seen as a quasi-harmonic series 
of oscillations of decreasing amplitude in the power spec- 
trum at wave iiumbers O.OlftMpc"^ < k < 0.4hMpc~^ 
l|Sugivamalll"995l : lEisenstein fc Hull 19981119991 1. In the two- 
point correlation function, the Fourier transform of the 
power spectrum, the B AO are visible a s a unique broad 
and quasi-gaussian peak (|Matsubajall2004l l . The BAO where 
flrst detected in the correlation function of the Luminous 
Red Galaxies (LRG) sample drawn from t he Slo an Digi- 
tal Sky Survey (SDSS) by Eiscnstci n et al.l (|2005h and the 
power spectrum of t he two-degree Fie ld G alajcy Redshift 
Surv e y (2dFGRS) bvlCole et al.l (|2005 ll fsee I Percival et all 



20071: ICabre fc Gaztanagal l2009l: J Gaztafiaga et al.1 12001 



20101 : 



Sanchez et al.' '20091: iPercival et ah! [2009i : iReid et al.l 
Kazin ct al. 2009, for more recent analyses). The acoustic 
scale inferred from the galaxy clustering is related to the 
sound horizon scale at the drag epoch, i.e. when the baryons 
where released from the photons. Because of the very high 
photon to baryon ratio, this ep och happe ned slightly later 
that the decoupling epoch (K omatsu et a l. 2009). As this 
scale depends solely on the plasma physics after the big bang 
and can be calibrated using CMB data, it is possible to use 
the BAO scale as a standard ruler. Measuring the apparent 
size of the BAO in the directions parallel and perpendicular 
to the line of sight, it is possible to measure the redshift 
dependence of the Hubble parameter H and the angular di- 
ameter distance Da and thu s constrain cosmologica l param- 
eters like wpE JBlake fc G lazcbrook 2003; Hu fc HaimanI 
2003:'Li nder et al.! '2003": Sco & Eiscnstcin 2003; Wan'g 2006; 
Guzik. Bernstein fc Smith .2007.: .Sep fc Eisenstein .2007.: 



Seo et al.l l2008l: IShoii, Jeong fc Komatsul l2009l: ISeo et al.1 



20091 ). 



Robust percent level constraints on iude re- 
quire measurement of th e BA O scale at sub -p ercent 
level (iNishimichi et ai] 120071 : lAngulo et all l2008l : 
IShoii. Jeong fc Komatsul [2009|). In order to realize the 
full potential of future BAO measurements as a cosmolog- 
ical probe, a detailed analysis of the possible systematic 
effects introduced by the non-linear evolution of den- 
sity fluctuations, redshift space distortions and bias is 
required. Recent studies based on the analysis of large 
volume N-Body simulations and on novel approaches 
to pert urbation theory (e.g. renorm alized perturbation 
theory, ICrocce fc Scoccimarroll2006al lal. have shown that 
non-linearities introduce shifts of up to few percent in the 
position of the BAO peaks with respect to the prediction 
of linear theory, both in the correlation function and in 
the power spectrum JSmith. Scoccimarro fc ShethI l2007l : 
ICrocce fc Scoccimarrd l2008l : ISanchez et al.ll2008l '). Redshift^- 
space distortions and bias might introduce additional shifts. 
When not properly taken into account, these distortions 
introduce systematic errors in the obtained constraints that 
can be larger than the expected measurement errors. This 
highlights the importance of using a model able to take in 
account these distortions with the accuracy demanded by 
the future surveys. 

The modelling of non-linear distortions in the BAO 
signal in the power spe ctrum has been t h e focus of 
much theoretical work JAngulo et all l2005l: iHuff et al 
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Crocce fc Scoccimarrd l2008l : ISeo et al.l l2008l : ISeo et al.1 
2009'). Most of these analyses have focused on modelling 
the BAO signal flltered out from the broad band shape of 
the power spectrum. This approach attempts to produce a 
purely geometrical test based on the BAO signal, but it has 
the disadvantage of discarding useful informa tion contained 
in the full shape of the power spectrum. I Sanchez et al.l 
((2003) showed that the correlation function is much less 
affected by scale dependent effects than the power spectrum 
and that its full shape can be described by a s imple model, 
originally proposed by ICrocce fc Scoccimarrd (|200S' ) , with 
the accuracy required for surveys of up to two orders of 
magni tude larger than present day samples. ISanchez et al.l 
(|2009l ) applied this model to the LRG co rrelation function 
measured by IC abre fc Gaztaiiagal ([20091) and performed 
a detailed analysis of the constraints on cosmological 
parameters, obtaining an improvement of roughly a factor 
two in the accuracy of the constraints on the dark energy 
equation of state with respect to the ones recovered by 
using only the BAO signal in the power spectrum. 

In this work we present a new model for the full shape 
of the power spectrum which is based on the mode l for 
the correlation funct i on oflCrocce fc Scoccimarrd (|2008r ) and 
ISanchez et al.l (|2008l . |2009| ) . We compare this model against 
the real and redshift space power spectra of the dark mat- 
ter and halo distributions obtained from an ensemble of 
large volum e, moderate resolut i on N-body simulatio ns (L- 
BASICC II lAngulo et all I2OO8I : ISanchez et ahl 120081 ). The 
outline of the paper is as follows: in section [5] we describe 
the set of simulations used in our work and provide a few 
technical details about the computation of the power spec- 
trum. In section [3] we introduce our theoretical model of the 
full shape of the power spectrum. In Section |4] we determine 
the range of scales in which this model is able to accurately 
describe the results drawn from the numerical simulations 
and recover unbiased constraints on the dark energy equa- 
tion of state. Finally section[5]presents our main conclusions. 



2 N-BODY SIMULATIONS AND THE 
COMPUTATION OF THE POWER SPECTRUM 

In this section we briefly describe the L-BASICC II e nsemble 
of N- Body simulations used in our analysis (see Ang ulo et al.l 
|2008| . for more details) and give some technical details about 
the methodology implemented to compute the power spectra 
from these simulations. A more detailed description can be 
found in Appendix 1X1 



2.1 The L-BASICC II N-Body simulations 

We use an ensemble of 50 moderate resolution, very large 
vol ume dark matter N-Body simulations c alled L-BASICC 
II (|Angulo et al.ll2008l : ISanchez et al.ll2008l ). These are ana- 
logues of the L-BASICC simulations used in lAngulo et al.l 
(2008) and represent the evolution of the dark matter den- 
sity fleld in a universe characterised with a flat ACDM 
cosmology, consistent with the constraints on cosmologi- 
cal parameters obt ained from the c ombin atio n of CMB and 
LSS in formation of ISanchez et al.l (|2006l ) and lSpergel et al.l 
(2007). The values of the cosmological parameters and other 
speciflcations of the simulations are listed in Table [l] 
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Table 1. Cosmological parameters and specifications of tfic 
L-BASICC II simulations. 



matter density 


f^m 


0.237 


baryonic density 


Hb 


0.041 


scalar spectral index 


"s 


0.954 


Amplitude of density 
fluctuations 


T8 


0.77 


Hubble constant 


Ho 


73.5kms-iMpc-i 


Number of particles 


TVp 


448^ 


Particle mass 


Mp 


1.75 X 10^2 h-iM© 


Softening length 


e 


200 h-^ kpc 


Comoving box side 




1340 h-^ Mpc 


Comoving volume 




2.41 h-^Gpc^ 


Total volume of the ensemble 




120 h-^Gpc^ 



Table 2. Number of haloes (N^^) and shot noise term (l/n) 
for the total halo sample and the three mass bins defined in 
section [JT] for redshift 0, 0.5 and 1. 





z = 


:0 


z = 


0.5 


z = 


: 1 


bin 


A^h 


1/n^ 


A^h 


l/n 


A^h 


l/n 


tot 


465903 


5164 


294204 


8178 


143531 


16764 


1 


262232 


9175 


151976 


15832 


69457 


34642 


2 


101825 


23630 


71551 


33628 


36852 


65291 


3 


101846 


23625 


70678 


34043 


37221 


64644 



" units of h-^ Mpc^ 

The position and the velocity of all the particles in 
the simulations were stored in three snapshots at redshifts 
z = 0, 0.5 and 1. Halo catalogues were constructed from the 
dark matter distributions at each redshift u sing a Friend-of- 
Friends (FoF) algorithm l|Davis et al-lllOSSl l. with a linking 
length parameter b = 0.2 and selecting all the haloes with 
more than 10 particles, which corresponds to a minimum 
halo mass of 1.75 x lO" h'^MQ. 

From each halo catalogue we extract three sub samples 
selected according to the following mass limits (in units of 
10^^/i-^Mq): 

1) M < 3.5, 

2) 3.5 < M < 5.95, 

3) M ^ 5.95 

at z = (selected as in lSanchez et al.ll2008h . 

1) M < 2.9, 

2) 2.9 < M < 4.65, 

3) M ^ 4.65 

at z = 0.5 and 

1) M < 2.6, 

2) 2.6 ^ M < 3.8, 

3) M ^ 3.8 

at z — 1. These limits were chosen in order to include about 
half of the total number of haloes in mass range 1 and the 
remaining equally divided between samples 2 and 3. The 
number of haloes (A^h) and the shot noise term (l/n), a scale 
independent poisson term arising from the discretization of 



the density field (see also Appendix [XJ , for each mass bin 
at each redshift are given in Table [2] 



2.2 Power spectrum computation and shot noise 

In order to compute the power spectra, the dark mat- 
ter particles or haloes in the simulations were assigned to 
a grid of 1008^ cells using the Triangular Shaped Cloud 
(TSC) as Mass Assignment Scheme (MAS). We then com- 
puted the Fourier transform of the obtained density field 
by a fast Fourier transform (FFT) algorithm using the free 
software FFTWj (Fas test Fourier Transform in the West, 
iFrigo fc Johnsonll2005l ). We then correct the amplitude of 
the Fourier modes for the effects of the MAS as in the 
first line of equation (|A5b[) . spherically average them in 
shells of thickness Afc = 2n/L = 0.0047 /iMpc"^ and sub- 
tract the shot noise contribution l/n. For this configura- 
tion the Nyquist wavenumber is kjq = 2.36/iMpc~^ and 
the computed power spectrum is exact for k < 67%feN = 
1.58/iMpc~^ (see appendix 1X1 for more details). To obtain 
the redshift space power spectrum we computed the appar- 
ent position of the dark matter particle or haloes converting 
their velocities along a single axis of the simulation into a 
displacement in comoving coordinates. 

We computed the real and redshift-space power spectra 
of the dark matter distribution (Pdm) at redshift 0, 0.5 and 
1, and the corresponding power spectra of the total halo 
catalogue and the three mass sub-samples, which we label 
as Ptot, Pii, P22 and P33 respectively. We also computed the 
cross power spectra for the three possible combinations of 
the mass bins at each redshift (P12, P13 and P23). 

The FoF algorithm, used to create the halo catalogues, 
is intrinsically exclusive: two haloes must be separated by 
a distance larger that the sum of their radii or they would 
be identified as a single more massive halo. This introduces 
an exclusion effect in the halo catalogues which is visible in 
both the correlation functions and the power spectra of these 
samples. In the correlation function the exclusion effect is 
clearly visible at small scales, where it becomes negative. As 
an example. Figure [l] shows the mean correlation function 
obtained from the total halo catalogue at 2 = in our en- 
semble of simulations: after reaching its maximum value at 
r ^ 2h~^ Mpc, it decreases converging to ^ = — 1 as r — >■ 0. 
Figure [2] shows the mean power spectrum with and without 
the shot noise subtracted (solid and long-dashed lines, re- 
spectively) for the same sample and redshift. The dot-dashed 
lines show their corresponding 1 — a variance. The shot noise 
amplitude is indicated by the horizontal dashed line. It is 
important to notice that when the shot noise is subtracted, 
the power spectrum becomes negative for k > l/iMpc~^. 
We have tested that this feature is independent on the shot 
noise amplitude or on the dimension of the grid used for the 
FFT. This clearly points to the fact that in the presence of 
the exclusion effect the noise is no more Poisson and possibly 
scale-dependent . 

The problem of the existence and impact of the ex- 
clusion effect in the power sp ectrum has been already ad- 
dresse d in previous analyses. ISmith. Scoccimarro fc ShethI 
(|2007|) point out that the exclusion effect may give rise to a 
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10 
X [Mpc/h] 



100 



Figure 1. Log-linear scaling of the mean real-space correlation 
function S,{x) (solid line) of the total halo sample from the ensem- 
ble of simulations. The dotted lines indicate the variance from the 
different realisations. The signature of the exclusion effect in the 
halo sample can be clearly seen at r < 2 h~^ Mpc. 



scale dependent noise term which may lead t o a misinterpre- 
tation of t he shape of the power s pectru m; ICasas-Mirandal 
l|2002h and iManera fc Gaztanagal (|2009h found that biased 
objects present non Poisson noise. Since the understanding 
and modelhng of the influence of the exclusion effect on the 
power spectrum is beyond the scope of this paper, we will 
not address specifically this problem here. Nonetheless, as 
will be explained section [S] our model of the shape of the 
power spectrum contains few free parameters that can par- 
tially absorb possible deviations from white noise. We will 
come back on this issue in section [S] 



3 MODELLING THE FULL SHAPE OF P(K) 

3.1 The model 

The increasing volume of the new large galaxy redshift sur- 
veys requires accurate models of the LSS observations in or- 
der to extract the maximum amount of information from the 
data without introducing systematic effects. The two stan- 
dard approaches used to model the power s pectrum and the 
corre lation function are the halo model fsee lCoorav &: ShethI 
12003 . for a r eview) , and the use of p erturbation theory tech- 
niques (see iBernardeau et al.l |2002| . for a review) . At red- 
shift 2 > 1, perturbation theory is able to provide a cor- 
rect description of the power spectrum of the dark mat- 
ter clustering at mildly non linear scales when terms of at 
least third order in the densi ty fluctuations are included 
l|jeong fc Komatsrj|2006l . l2009l ). However, none of these ap- 
proaches is enough to model the full shape of the power 
spectrum in the non-linear regime or at smaller redshifts. 

In recent years, there has been substantial progress re- 
garding the theoretical understanding of non-linear evolu- 




0.01 0.10 1.00 

k [h/Mpc] 

Figure 2. Mean power spectrum with and without the shot noise 
subtracted (solid and long-dashed lines, respectively) from the en- 
semble of simulations and their corresponding l-cr variance (dot- 
dashed lines) for the total halo sample. The amplitude of the shot 
noise is indicated by the horizontal dashed line. To enhance the 
negative part of the power spectrum at small scales, we use linear 
vertical axis for P{k) < 200hMpc~^ instead of logarithmic. 



tion using perturbation theory. Have been proposed sev- 
eral methods which are based on a partial resummation of 
high order terms both in Eulerian space as renormaliz ed 
perturbation theory (RPT, Crocce fc Scoccimarrdl2006al lbl: 



IBernardeau. Crocce fc Scoccimarrdl2008l ) or in Lagrangian 
space (Lagrangian Perturbation theory: LPT, iMatsubaral 
i2008aH bl). Similar approaches are based on the renor- 
maliz ation group equations (jMatarrese fc Pietronil 120071 , 
l200a ) or on the renormalization of the parameters of the 
model in order match observables (IMc Donald 200f|, 120071 : 
ISmith. Hernandez- Monteaeudo fc Seliakl l2009l ) . A different 
approach to PT is based on breaking the infinite hierar- 
chy of n-point statistics by introducing p hysically moti- 
vated closure equations at the desired orde r |Pietronill2008l : 
Taruya & Hiramatsu 2008; Taruya et al. 20^01^ 

Within the theoretical framework of RPT, the non- 
linear power spectrum of the density fluctuations can be 
computed as the sum of two terms containing different 
physical information. The flrst of these contributions con- 
tains all the terms in the perturbation theory expansion 
of P(k, z) that are proportional to the initial linear theory 
power spectrum Ph{k) in the so called renormalized propa- 
gator G{k, z), which represents a non-linear correction to the 
growth factor. The second contribution to P{k, z) groups all 
the remaining terms into the mode coupling power spectrum, 
Pmcik, z), that contains the power arising from the coupling 
of different Fourie r modes. The result ing power spectrum 
can be written as IjCrocce fc Scoccimar ro 2006a b) 



P{k, z) = G(fc, zfP^(k, z) -f PMc(fc, z), 



(1) 
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where the redshift evolution of Ph{k, z) is given by the linear 
theory growth factor D(z). 

In the large-fc limit the propagator is accurately de- 
scribed by a Gaussian damping 



G{k, z) = exp 



f k 



\V2k, 



where the damping scale fc* is given by 
1 



k^z) 



67r2 



dfcPL(A:, z] 



-1/2 



(2) 



(3) 



Equation © also provides a good aproximation for the prop- 
agator for small values of k. 

The mode coupling power spectrum is given by the sum 
of an infinite number of terms Pn Mc{k) ordered by the num- 
ber n of initial modes coupled. The contribution of all the 
terms combining n modes becomes important at progres- 
sively smaller scales as n increases. The first term is thus 
P2 MC (fc), the power generated by the coupling of two modes, 
which can be approximated by the standard PT term 



-Plloop(fc) - -^ 



d3qlF2(k-q,q)l2pL(|k-qj)PL(g), (4) 



where P^Ck, q) is the second order kern el of perturbation 
theory (eq. 45 in iBernardeau et al.l 120021 ). In the literature 
the power in equation Q is usually reffered to as P22{k). 

Panel a of Figure [3] shows the linear theory power spec- 
trum (solid line) and Piioop(fc) (dashed line) for the cos- 
mological model of our ensemble of simulations. Panel b 
shows Phj k) divided by a referen ce power spectrum with- 
out BAO (JEisenstein fc Hu" 1998) and Piioop(fc) divided by 
the corresponding term derived applying equation Q to the 
smooth power spectrum. Piioop(fc) also contains oscillations, 
although with smaller amplitude and, more importantly, out 
of phase with respect to the ones in Pi^{k). When these two 
terms are summed as in equation ([l}, the BAO are shifted 
towards smaller s cales with respect to the one s in the linear 
p ow er spectrum |Crocce fc Scoccirnarrd uOOSl 'l . 

ICrocce fc Scoccimarrd (2006b|) showed that the pre- 
dictions from RPT give an accurate description of the 
real space dark matter power spectrum measured from N- 
body simulations without the need to tune a free param- 
eter. Despite this success, a drawback of this formalism 
is that there is no straightforward way to include the ef- 
fects of redshift space distortions and bias, making its direct 
application to observa t ional datasets impossible. However, 
ICrocce fc Scoccimarrd (|2008l ) proposed a model for the large 
scale correlation function motivated by the RPT formalism. 
In this ansatz the correlation function is given by 



CNL(r) = 6^UL(r) 



-{ktry 



+ AMC^i{r)^i'\r: 



(5) 



where b, fe* and Amc are treated as free parameters, and the 
symbol ® denotes a convolution. Here ^l is the derivative of 
the linear correlation function and ^^ (^) is defined by 



e^w 



-^c.w - ^ 



Pi^{k)ji(kr)kdk, (6) 



with ji{y) denoting the spherical Bessel function of order 
one. The second term in equation © corresponds to the 
leading order contribution to ^mc from the one-loop ap- 
proximation of the mode coupling power of equation (|4|. 




0.01 



0.1 
k [h/Mpc] 



Figure 3. Panel a: Linear theory (solid line) and first loop 
(dashed hne) power spectra. Panel b: ratio between the power 
spectra of panel a and a reference power spectrum without oscil- 
lations. The first loop contribution PnoopC^) shows small oscilla- 
tions out of phase with respect to Piin{k), generating a net shift 
of the BAO peaks when summed. 



I Sanchez et al.l (|2008l ) compared this model against the re- 
sults of N-body simulations and found that it is able to 
give an accurate description of the full shape of the cor- 
relation function, including the effects of bias and redshift 
space distortions, for volumes up to two orders of magni - 
tude larger than present day datasets. I Sanchez et al.l (|2009l ) 
successfully used this model to obtain constraints on cosmo- 
logical parameters from the correlation function of a sample 
of luminous red galaxies (LRG) drawn from the data release 
6 (DR 6) of the SDSS as measured bv lCabre fc Gaztafiagal 
(|2009l ). 

In our analysis, we follow the same approach and model 
the non-linear power spectrum as 

P(fc, z) = b^ (e"<'/'=*''Pii„(fc, z) + AMcPuoopik, ^)) , (7) 

and treat b, fc* and ^mc as free parameters. In the next 
section we will show that this model allows to obtain un- 
biased constraints on the dark energy equation of state 
parameter by accurately describing the full shape of the 
power spectrum measured in real and r edshift space. Re- 
cently [Croccei^Cabr£^^£azt^ag3 (|201G| ) pointed out that, 
according to equation ^ , fc* scales with the inverse of the 
growth factor D{z)~^ . 

Piioop(fc) does not have the same shape as P2Mc(fc): 
the latter in fact decreases faster than the former and at 
k ~ 0.15 — 0.2/iMpc~^ it is roughly 1.5 — 2 times smaller. 
But at those scales the amplitude of the three mode cou- 
pling, P3Mc(fc), is already around 1/4 — 1/2 of P2Mc(fc) 
and shoul d be included in the model (se e Figure [3] here and 
figure 1 in lCrocce fc Scoccimarrdl2006al . for a comparison). 
Piioop(fc) is thus somewhat larger than P2Mc(fc) + P3Mc(fc) 
and the difference becomes more and more important with 
increasing wave number. We will come back to this again in 
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section |4l where we will discuss the range of scales in which 
the model of equation Q can be applied to a measurement 
of the power spectrum. 

We expect our model to be less efficient at describing the 
shape of the power spectrum in redshift-space than in real- 
space. In fact, we do not include redshift space distortion 
explicitly and we let the free parameters compensate some 
of its effects. Since the scale dependence of the redshift space 
distortions is stronger for the dark matt er than for the haloes 
ljScoccimarrdl2004l : lAngulo et al.ll2008l ). the lack of a model 
for them will be particularly visible in the former case. 



4 MODEL IN PRACTICE AND DISCUSSION 

In this section we test whether or not the model described 
in section [5TT] can give an accurate description of the power 
spectrum including the effects of non-linear evolution, red- 
shift space distortions and bias. In section 14.11 we describe 
the test we implement to determine if the model returns un- 
biased constraints on the dark energy equation of state pa- 
rameter iude. Sections 14.2)14.41 describe the results obtained 
when each of these scale dependent effects is included in the 
measurement of P{k). 



4.1 Testing the model 

We now test whether or not the model for the power spec- 
trum described in section lSTTI returns unbiased constraints on 
the dark energy equation of s tate parameter mpE. For this, 
we foUow lAngulo et al.l (|2008l ') and lSanchez et al.l l|2008l ') and 
consider the simple case in which we assume the values of 
all other cosmological parameters to be known, and analyze 
the constraints on wde only. 

When measuring the power spectrum from an observa- 
tional dataset, it is necessary to assume a fiducial cosmology 
in order to map the observed galaxy redshifts and angular 
positions into comoving distances. This choice has an im- 
pact on the results obtained. The use of a value of wde dif- 
ferent from its true value modifies the measured separations 
between the objects in the sample, changing the shape of 
the measured power spectrum. However, for small deviations 
away from the true equation of state, the alteration in the 
measured power spectrum can be represented by a rescaling 
of the wavenumber from fctruc to fcapp. This cha nge can be 
encap sulated in a 'stretch' factor q, defined by (JHuff et al.l 



encapi 
I2OO7I ') 



fcti 



(8) 



We analyse the constraints on the stretch parameter 
obtained from the mean power spectra of the L-BASICC 
II simulati ons, which can be tr anslated into constraints on 
«;de as in lAngulo et al.l (|2008r ) . All quoted allowed ranges 
for the constrained parameters correspond to the 68% con- 
fidence level according to the variance from the ensemble of 
simulations. 

We explore the parameter space defined by 
9 — (ki,, 61, Amc, ol) us ing a Markov chain Monte Carlo 
(MCMC) te chnique JGilks. Richardson fc Spiegelhalteij 
1 19961 : IChristensen fc MeveJ I2OO0I '). We assume that the 



likelihood function follows a Gaussian form 

1 



C oc exp 



-yio) 



where 



x\o) 



(d 



t{e)fc-\A- 



t{e)) 



(9) 



(10) 



is the standard x^, in which d is an array containing the 
fitted data, in our case the computed power spectrum, t{9) 
contains the model computed for a given set of parameters 
9 and C is the covariance matrix of the measurement. 

We computed the covariance matrix C for the dark mat- 
ter and the halo power spectra from the L-BASICC II sim- 
ulations as 



JVroal 






-y E(^'(^')-nfcO)(^<fe)-^fe)), (11) 



where Piiki) corresponds to the measurement of the power 
spectrum at the i-th fc-bin in the Z-th realisation and P{ki) 
corresponds to the mean power spectrum from the ensemble 
at the same wavenumber. 

Figures U and [S] show the correlation matrices 
CijI \jCiiCjj of the power spectra of the real-space dark 
matter and total halo samples respectively. In both figures, 
the three panels correspond, from left to right, to z = 0, 
0.5 and 1. In all cases we observe that the correlation be- 
tween different modes, due to non linear mode coupling, is 
stronger at 2: = and decreases at increasing redshift. At 
2 = the correlation matrix of the total halo sample shows 
strong correlations for fc > 0.1/iMpc~^. In order to obtain 
robust constraints on a these correlations must be included 
when fitting our model to the results from the L-BASICC 
simulations. For a more detailed studies on the covariance 
of the power spectrum using theoreti cal model and large 
numerical simulations sec Hamilton, Rimes fc Scoccimarrd 
^20061 ): ISmithI (|2009l ): lTakahashil (|2009l ). 

On our parameter space we implemented flat priors 
given by 

• 0/iMpc"V fc* < 0.35/iMpc"\ 

• s; Amc < 10, 

• 0.5 s; Q < 1.5. 

We analytically marginalize the bias parameter h over an in - 
finite flat prior using equation F2 in lLewis fc Bridle! l|2002h . 
In order to obtain an estimate of the amplitude of the model 
power spectrum, we compute a value of h and its variance 
by maximizing the likelihood function of equation Q while 
the other parameters are kept fixed to their mean values 
obtained from the MCMC. 



4.2 Non-linear evolution 

The circles in Figure|n]show the mean real-space dark matter 
power spectra from the L-BASICC II simulations at z = 0, 
0.5 and 1 (panels a, h and c respectively). Dot-dashed lines 
show the corresponding variances from the estimates in the 
different realizations. In order to highlight the signature of 
the acoustic oscillations, Figure [7] shows, with the same sym- 
bols, the power spectra divided by a smooth linear theory 
power spect rum without BAO com puted using the fitting 
formulas of lEisenstein fc Hul (|l998i V The comparison with 
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Figure 4. Correlation matrix Cij/ ^CaCjj, where dj are the elements of the covariance matrix C, for the real-space dark matter power 
spectra at 2: = (left), z = 0.5 (center) and 2 = 1 (right). At z = 1 the correlation matrix is close to diagonal, but at lower redshifts, when 
non-linearities become increasingly important, a significant correlation between different modes appears and becomes more important 
for larger scales, i.e. smaller k, as the redshift decreases. A similar result is obtained for the redshift-space power spectrum. 



haloes z = 



z = 0.5 



z=1.0 



u 

CL 
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Figure 5. Same as in Figure |4] but for the total halo samples at 2; = (left), z = 0.5 (center) and 2 = 1 (right) in real space. As for 
dark matter, the mode coupling becomes more important at larger scales as the redshift decreases. A similar result is obtained for the 
redshift space catalogue and for the other halo catalogues. 



u 
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Figure 6. Mean power spectra computed from the simulations (circles for real space and triangles for redshift space), their variance (dash- 
dotted lines) and the model power spectrum as obtained through the fitting (solid lines) as function of the wavenumber in comoving units 
for the dark matter catalogue at redshift 0, 0.5 and 1.0 (left, centre and right respectively), in log-log scaling. The maximum wavenumber 
used for the fit is fcmax = 0.15/iMpc~^ and is indicated by the vertical arrow. 
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Figure 7. Power spectra of Figure [6] divided by a smooth reference power spectrum llEiaenstein &: Hulll998f) in order to enhance the 
BAO oscillations, in linear scahng. For comparison the linear power spectrum, divided by the smooth power spectrum, is shown (dashed 
hnes). The upper horizontal hnes are drawn to the values corresponding to the Kaiser boost factor (section 14.31 and Table |4]|. The smooth 
power spectrum is the same for all the plots. The maximum wavenumber used for the fit is fcmax = 0.15 h Mpc"'^ and is indicated by the 
vertical arrow. 



the linear power spectrum (dashed hnes) shows that non lin- 
earities change the broad band shape of the power spectrum 
al small scales and damp the BAO feature. We explored the 
constraints on the stretch parameter obtained by applying 
the model described in section [HTT] to these measurements. 

The hnes in the panels a, h and c of Figure [S] show the 
mean value of a obtained from the measurements at z = 0, 
0.5 and 1 (solid hnes), together with its correspondent 68% 
and 95% confidence levels (dot-dashed hnes), as function of 
fcmax, the maximum value of k included in the analysis. As 
smaller scales (larger values of fcmax) are taken into account, 
the width of the allowed region for a decreases due to the 
larger number of modes included in the fit. At z = 0, the 
mean value of a remains consistent with 1 at a 1-cr level 
for fcmax < 0.16 /iMpc-\ with a = 0.999 ± 0.007. At higher 
redshifts the value of fcmax for which this holds increases, 
due to the smaller impact of non-linearities. The first row 
of the upper part of Table |3] lists the obtained values of a 
for fcmax = 0.15 /iMpc~^. The solid hnes in Figures [6] and [7] 
show the model power spectrum of equation computed 
using the mean values of the four parameters obtained for 
this range in fc. This model is able to accurately describe the 
effects of non-linear evolution in both the broad-band shape 
of the power spectrum and the damping of the acoustic os- 
cillations, up to the maximum value of fc included in the fit, 
which is indicated by a vertical arrow. 

The dashed contours in Figure [9] show the two- 
dimensional marginalized constraints in the a — fc* plane 
obtained for fcmax = 0.15 /iMpc~^ at 2 = 0, 0.5 and 1 
(panels a, b and c respectively). The contours correspond 
to Ax^ ~ 2.3 and 6.17 which, assuming a two-dimentional 
gaussian likelihood, are equivalent to the 68% and 95% confi- 
dence levels. While the constraints on the stretch parameter 
are very tight, there is a wide allowed region for fc* whose 
mean value is larger than the theoretical prediction of Equa- 
tion ([3]). The mean value fc* shows a tendency to increase 
with redshift, as the BAO feature is less damped, but the al- 
fowed range for this parameter is too large to compare this 
evolutions with D{z)~^ . The wide allowed range for fc* is 
caused by a strong degeneracy between this parameter and 



Table 3. Mean values of a and their I-ct c.l. as recovered 
from the mean dark matter and halo samples power spectra of 
the three redshifts outputs of the L-BASICC II simulations for 
kmax = 0.15/iMpc~^. The upper part is for real space, the lower 
for redshift space. See Section [2.2l for the definition of the different 
samples listed in the second column. 







z = 


2 = 0.5 


z = 1 




DM 


0.999 ± 0.007 


1.001 ±0.006 


1.001 ± 0.006 




tot 


1.001 ±0.011 


0.995 ±0.011 


0.994 ±0.014 




11 


0.999 ± 0.014 


1.006 ±0.014 


0.999 ±0.017 




22 


1.006 ±0.017 


1.006 ±0.013 


1.005 ± 0.023 


Olreal 


33 


1.003 ±0.016 


0.997 ±0.012 


0.987 ±0.016 




12 


1.002 ±0.013 


1.004 ±0.01 


0.987 ±0.015 




13 


1.002 ±0.01 


0.991 ±0.012 


0.991 ±0.018 




23 


0.995 ±0.013 


1.004 ±0.014 


1.005 ± 0.013 




DM 


0.994 ± 0.007 


1.001 ±0.007 


1.005 ± 0.008 




tot 


1.002 ±0.012 


0.997 ±0.01 


0.987 ±0.012 




11 


0.999 ±0.016 


1.008 ±0.016 


0.998 ± 0.016 




22 


1.007 ±0.02 


1.01 ±0.018 


1.006 ± 0.026 


"red 


33 


0.998 ±0.017 


0.994 ±0.014 


0.981 ±0.015 




12 


1.002 ±0.014 


1.006 ±0.012 


0.982 ± 0.014 




13 


1.008 ±0.011 


0.993 ±0.012 


0.978 ± 0.018 




23 


0.997 ±0.015 


1.006 ±0.014 


1.006 ± 0.014 



Amc- This can be seen in Figure [TOl which shows the two- 
dimensional constraints in the Amc — fc* plane for the dark 
matter power spectrum at 2: = 0, 0.5 and 1. The degeneracy 
arises because it is possible to provide a good description of 
the overall shape of the power spectrum by compensating 
an increase in the damping of the first term of equation ((Tjl 
(a decrease of fc*) by increasing the amplitude of the mode 
coupling contribution (using a higher value of Amc)- Be- 
sides this, the value of ^mc obtained at 2; = is about 30% 
smaller than the expected, confirming that the Piioop(fc,2) 
is somewhat bigger than the sum of the RPT two and three 
mode coupling. This makes the value of fc* slightly larger 
than the theoretical value of equation ^. In Figure [TU] it is 
also noticeable that the values of Amc increase with redshift 
due to the smaller relative amplitude of Piioop(fc), which de- 
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Figure 8. Mean value of a obtained from the MCMC (blue solid line for real space, red long-dashed line for redshift space), 1, 2 — ct 
confidence level (dotted lines) as function of fcmax for dark matter (real space: panels a, b and c; redshift space: panels d, e and /) and 
the total halo catalogue (real space: panels g, h and i; redshift space: panels j, k and /) at z=0, 0.5, 1 from left to right. The horizontal 
dashed lines indicate o = 1. 
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Figure 9. 1 and 2 a contours of likelihood map in the o — fc* plane obtained from the MCMC for dark matter at redshift (left), 0.5 
(centre), 1 (right), kmax = 0.15/iMpc~^. Background areas within dashed lines are for real space, foreground ones within solid lines are 
for redshift space. 
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Figure 10. 1 and 2 a contours of likelihood map in the A^c ~ ^* plane obtained from the MCMC for dark matter at z=0 for 
kmax = 0.15/iMpc~^. Background areas within dashed lines are for real space, foreground ones within solid lines are for redshift space. 




creases as the growth factor squared, with respect to the 
linear power spectrum. 



4.3 Redshift space distortions 

In galaxy surveys distances are inferred from the measured 
redshifts. The observed redshift of a galaxy is given by the 
sum of the contribution from the cosmological expansion and 
the Doppler shift due to its peculiar motion along the line 
of sight. This gives rise to a difference between the real and 
the apparent position of a galaxy, leading to an increase in 
the amplitude and a change of shape in the measured power 
spectrum Ps{k) with respect to its real-space counterpart. 

On very large scales, peculiar velocities are dominated 
by the coherent flow of matter towards over-dense regions. 
In linear theory, under the plane parallel approximation, this 
produces a scale-independent boost in the amplitude of the 
large scale power spectrum given by (JKaiseill 19871 ) 






(12) 



where /3 = f/b, b is the bias factor and / is the logarithmic 
derivative of the growth factor D with respect to the scale 
factor a 



/ = 



dlnD 
din a 



(13) 



The linear theory description of equation (1121) is only 
valid asymptotically on extremely large scales CScoccimarrd 
l2004l : lAneulo et aDl2008l : I Jennings. Baueh fc Pascolill2010l l. 
On smaller scales, peculiar velocities are dominated by the 
random motions inside virialized structures which make 
bound haloes to appear elongated along the line of sight 
when mapped in redshift-space, an effect commonly known 
as "fingers of god" . This causes a damping of the power spec- 
trum that introduces deviations from the simple description 
of equation p2l) even at scales k > 0.03 hMpc~^ . 

Figure |6] shows the mean redshift-space dark matter 
power spectra (triangles) from the L-BASICC II simulations 
at z = 0, 0.5 and 1 (panels a, b and c respectively), together 
with their correspondent variances from our ensemble of sim- 
ulations (dot-dashed lines). Figure [7] shows the same power 



spectra divided by a smooth linear theory power spectrum 
without acoustic oscillation s computed using the fitting for- 
mulas of lEisenstein fc Hul (|l998ll . Since the smooth power 
spectrum is the same both in real and redshift-space, it is 
possible to see both the increase in amplitude and the change 
in shape towards smaller scales in the latter case. In this 
section we test if the parameterization of the model of equa- 
tion (O contains enough freedom to take into account these 
distortions. As stated in section [3Tl1 we can anticipate that 
the the performance of the model will be worse in this case 
than when dealing with real-space measurements, as we do 
not include explicitly the effect of redshift space distortions. 

The dashed lines in the panels d, e and / of Figure [8] 
show the mean value of the stretch parameter obtained by 
applying the model of equation (O to the mean redshift- 
space power spectra as a function of fcmax at z = 0, 0.5 and 
1, respectively. The dot-dashed lines indicate the correspon- 
dent 68% and 95% confidence levels. At z = 1, also when 
dealing with redshift-space information the model is able to 
recover constraints on a consistent with one at the 1 — ct level 
for k < 0.18 hMpc~^ . The constraints degrade at lower red- 
shifts as non-linear redshift space distortions become more 
important and, at 2; = 0, our results are only marginally con- 
sistent with a = 1. The first row of the lower part of Table [3] 
lists the constraints on a obtained for fcmax = 0.15/iMpc~^. 
In all cases, the allowed region for this parameter increases 
with respect to the real-space case. In all the panels of Fig- 
ures [6] and [7] the solid upper (orange) line shows the model 
power spectrum in redshift space with the parameters fixed 
to the mean value obtained from the MCMC. 

As described in section 14.11 we compute the value of b 
and its variance fixing the other parameters to their mean 
values as obtained from the MCMC. Since this represents 
the large scale linear bias, we can calculate the Kaiser boost 
factor simply as: 



Sflt 



bj_ 

b2' 



(14) 



with b and bs the real and redshift-space linear bias. The 
values of Sat obtained for the dark matter case are listed, 
together with the theoretical values Sun, in Table [l] The 
three values of 5'iin are also indicated in Figure [7] with the 
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Table 4. Effective linear bias (bcft) computed from the halo model 
prescription, values of b obtained from the fit in real space as 
described at the end of section [4JT] theoretical linear Kaiser boost 
factors Siin (equation [T2{ and Kaiser boost factors obtained from 
the fit {Sfit equation I14|l at all redshift, both for dark matter 
and the halo catalogues. All the ranges correspond to the 68% 
confidence level. 





Ka 


b 


•S'lin 


•Sfit 


DM 


1 


1.003 ±0.003 


1.34 


1.32 ±0.01 


tot 


1.89 


1.89 ±0.01 


1.17 


1.17 ±0.01 


11 


1.63 


1.65 ±0.01 


1.20 


1.20 ±0.02 


22 
^ = ° 33 


1.90 


1.84 ±0.01 


1.17 


1.17 ±0.02 


2.54 


2.49 ±0.01 


1.12 


1.12 ±0.02 


12 


1.76 


1.75 ±0.01 


1.18 


1.17 ±0.02 


13 


2.04 


2.04 ±0.01 


1.16 


1.17 ±0.01 


23 


2.20 


2.12 ±0.01 


1.14 


1.15 ±0.02 


DM 


1 


0.999 ± 0.003 


1.56 


1.55 ±0.01 


tot 


2.73 


2.65 ±0.01 


1.18 


1.18 ±0.01 


11 


2.31 


2.33 ±0.01 


1.22 


1.22 ±0.02 


— I 


2.59 
3.42 


2.57 ±0.02 
3.39 ±0.02 


1.19 

1.14 


1.19 ±0.02 
1.14 ±0.02 


12 


2.45 


2.45 ±0.01 


1.20 


1.21 ±0.02 


13 


2.81 


2.79 ±0.01 


1.18 


1.18 ±0.01 


23 


2.98 


2.95 ±0.02 


1.17 


1.17 ±0.02 


DM 


1 


0.996 ±0.003 


1.69 


1.70 ±0.01 


tot 


3.87 


3.80 ±0.02 


1.15 


1.16 ±0.01 


11 


3.27 


3.37 ±0.02 


1.18 


1.20 ±0.02 


1 22 
' = ^ 33 


3.58 


3.67 ±0.03 


1.17 


1.18 ±0.03 


4.62 


4.73 ±0.03 


1.13 


1.14 ±0.02 


12 


3.42 


3.50 ±0.02 


1.17 


1.18 ±0.02 


13 


3.89 


3.98 ±0.02 


1.15 


1.15 ±0.02 


23 


4.07 


4.18 ±0.03 


1.14 


1.15 ±0.02 



upper dotted lines. As can be seen the agreement between 
the theoretical and the recovered values are excellent at all 
the redshifts. 

The shaded contours within solid lines in Figure [9] 
show the two-dimensional 68% and 95% marginalized con- 
straints in the a — fc* plane obtained from the mean redshift- 
space power spectrum from the L-BASICC II simulations for 
fcmax ~ 0.15ftMpc~^ &t z — (panel a), 0,5 (6) and 1 (c). 
As in the case of real-space information, there is no degen- 
eracy between these parameters. Redshift space distortions 
increase the damping of the BAO signal. This is reflected 
in the mean values of fc* being systematically lower than 
the ones obtained from real-space data. Figure [TUl shows the 
correspondent constraints in the fc* — ^mc plane. The scale- 
dependent effects introduced by redshift-space distortions 
cause these two parameters to follow a different degeneracy 
than in the real-space case although the qualitative behavior 
is maintained. 



4.4 Halo bias 

In previous sections we compared the model of equation ((?]) 
against the clustering of the dark matter distribution. How- 
ever, when dealing with real observational data sets, only the 
clustering pattern of the galaxy distribution can be used to 
obtain cosmological information. The relationship between 
the galaxy and matter density fields can be highly non- 
trivial, adding an extra complication to the interpretation 



of clustering measurements in terms of constraints o n cos- 
mological parameters (e.g. see lSanchez fc Cole 1 12008). 

A general galaxy sample is composed of central and 
satellite galaxies. The former are assumed to lie at or near 
the centre of the dark matter haloes while the later can 
be associated with the sub-structure present in more mas- 
sive haloes that already contain a central galaxy. In red- 
shift space, the signature of the fingers-of-god effect can 
be associated with the contribution of the satellite galaxies, 
which then are responsible for the most imp ortant non-linear 
redshift space distortions. ,Tegmark et al.l (2004) analysed 
the SDSS DR2 com pressing the finge rs-of-god into isotropic 
structures. Recentlv lReid et al.l (|2009l ) proposed a method to 
reconstruct the underlying halo density field from a galaxy 
sample by identifying fingers-of-god like structures and re- 
placing them by one single halo. This procedure effectively 
eliminates the contribution of the satellite galaxies, which 
would correspond to the one-halo term in the halo-model 
formalism, from the measured power spectrum. This mini- 
mizes the impact of non-linear redshift-space distortions in 
the shape of the measured P{k). Hence, in order to extract 
cosmological information from these measurements, a model 
of the full shape of the halo power spectrum is required. 
iReid et all (|2009l ') also proposed a model, calibrated against 
numerical simulations, that relates the halo power spectrum 
to the linear theory prediction for the dark matter distribu- 
tion. In this section we compare the model of equation ((TJ 
against the halo power spectra of the L-BASICC II ensemble 
of simulations for different halo samples. 

According to linear theory, a simple relation is expected 
to hold between the dark matter and halo power spectra 
characterized by a scale-independent bias factor bi 



Phh(fc)=6?PDM(fc). 



(15) 



The results of numerical simulations have shown that 
this simple picture is only valid on very large scales 
and that, due to the effects of non-linear evolution, halo 
bias is a strong function of both scale and halo mas s 
(|Smith. Scoccimarro fc ShethI l2007l : lAngulo et all 120081 '). 
These distortions must be carefully modelled in order to ob- 
tain unbiased constraints on cosmological parameters from 
the full shape of the halo power spectrum. 

Figure [TT] shows the mean power spectra at z = of the 
total halo sample of the simulations, Ptot(fc) (panel b), the 
power spectra of the three mass bins described in section ETIl 
Pii(fc), P22{k) and P33(fc) (panels a, d and g), and their 
respective cross power spectra, Pi2(fc), Pi3(fc) and P2-i{k) 
(panels c, e and /). To increase the dynamical range of the 
plot, these power spectra have been divided by the same 
non-wiggle linear-theory power spectrum as in Figure[7| Cir- 
cles correspond to the real-space measurements while their 
redshift-space counterparts are shown by triangles, with the 
dot-dashed lines representing their respective variances from 
the ensemble of simulations. 

We obtained constraints on the parameter space defined 
in section 14.11 for different values of fcmax by applying the 
model of equation (O to our measurements of the halo auto 
and cross power spectra. As example. Figure [8] shows the 
constraints on a as a function of fcmax obtained from Ptot(fe) 
in real (panels g, h and i) and redshift-space (panels j, k and 
/) for z = 0, 0.5 and 1, respectively. The dot-dashed lines 
correspond to their 68% and 95% confidence levels. Due to 
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Figure 11. Mean power spectra computed from the simulations (circles for real space and triangles for redshift space), their variance 
(dash-dotted lines), the model power spectrum as obtained through the fitting (solid lines) and the linear power spectrum (dashed lines) 
as function of the wavenumber in comoving units at redshift for the total catalog (tot ), the three mass bins (11, 22 and 33) and the 
three cross mass bins {12, 13 and 23) divided by a smooth reference power spectrum iJEisenstein &: Hull 19981 ). Lower horizontal lines 
indicate the linear bias, b^ for the halo mass bin as computed from the halo catalogue prescription, the upper one indicate b^Si^^ (Table 
|4)|. The maximum wavenumber used for the fit is fcmax = 0.15/iMpc~ and is indicated by the vertical arrow. 



the lower number densities of the halo samples, the variances 
of these power spectra are larger than for the dark matter 
case and increase with increasing redshift. This leads to an 
increase in the allowed regions for a with respect to the 
ones obtained using Pumik). Despite this difference, these 
results show a similar qualitative behavior to the constraints 
obtained from the dark matter power spectrum, with a mean 
value consistent with a = 1 up to fc < 0.15/iMpc~^. The 
constraints on a obtained from the fit at different redshifts 
for fcmax = 0.15 /iMpc~^ are listed in TableO The solid lines 
in Figure [TT1 show the model power spectrum of equation ((7]) 
computed using the mean values of the four parameters ob- 
tained for this value of fcmax- 

The absence of a non-linear contribution from the 
fingers-of-god effect makes the scale dependence of red- 
shift space distortions much less significant for the 



halo power spectrum than for t h e dark matte r case 
IJSmith. Scoccimarro fc ShethI l2007l : JAngulo et all l2008l : 
iMatsubara l2008bl ) . This can be clearly seen in Figure ITTl 
where the differences in the shape of the real and redshift- 
space power spectra are small. This is reflected on the con- 
straints on the stretch parameters, which in the redshift- 
space case show a similar behavior to those obtained using 
the real-space measurements. 

Figure [12] shows the two-dimensional constraints in the 
a — fc* plane obtained using the mean auto and cross power 
spectra for the different halo samples at z = 0. The dashed 
lines show the 68% and 95% confidence limits obtained from 
real-space information while the solid lines correspond to the 
redshift-space case. The difference between the constraints 
obtained in the two cases is less significant than in the dark 
matter case shown in Figure (9] The constraints on fc* are 
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Figure 12. 1 and 2 a contours of the likelihood map obtained from the MCMC for the parameters a and fc* for all the haloes catalogues 
at redshift for fcmax = 0.15/iMpc~^. Background areas within dashed lines are for real space, foreground ones within solid lines are 
for redshift space. 



broad due to the degeneracy between this parameter and 
^Mc; nonetheless fc* seems to prefer smaller values in red- 
shift space, in accordance with the more important damping 
of the BAO features. 

In contrast to the behavior found in the analysis of the 
dark matter power spectrum, when applied to the different 
halo samples the ability of the model to obtain unbiased 
constraints on the stretch parameter degrades with increas- 
ing redshift and at z = 1 they are only marginally consistent 
with a = 1. This is related to the mass resolution of the L- 
BASICC II simulations: the smallest halo mass that can be 
resolved corresponds to Mhaio = 1-75 x IO^^/i'^Mq. While 
at 2 = these haloes are moderately biased tracers of the 
underlying dark matter distribution, at z = 1 they are much 
rarer objects that reside in very dense, and thus highly non- 
linear, regions. Hence, an accurate description of the shape 
of the power spectrum of such high mass objects at high 
redshifts requires a more detailed model of the non-linear 
distortions than that of equation ((7|. 



The shape of the cross power spectra Pi2(fc), Pis{k) 
and P2-i(k) can also be accurately described by the model 
of equation ^. As can be seen in Table |3] in all cases the 
constraints on a derived from these measurements are in 
good agreement with a = 1 and are tighter than the ones 
derived from their correspondent auto power spectra. 

As stated at in Section [2. 21 the halo catalogues are af- 
fected by the exclusion effect introduced by the Friend-of- 
Friends algorithm used to identify them. To account for a 
possible difference in the large scale amplitude of the shot- 
noise as could be caused by this effect, we have repeated the 
analysis adding a constant term to the model of equation 
which we allowed to vary over a wide flat prior. Our findings 
indicate that such a term is degenerate with fc* and AmCj 
while the constraints of a are not affected. 

We recover the values of the bias h for all the power spec- 
tra in Figure [TT] as indicated in section 14.11 The real space 
biases are listed in Table |4] and 6^ is shown by lower dotted 
lines in Figure 1111 As a comparison we also list the theo- 
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retical values of the effective bias bcff for the corresponding 
halo samples c omputed using the mass function proposed by 
Ijenkins et ali (2001 ) and the halo bias prescription obtained 
from the spherical collapse model (jSheth. Mo fc TormerJ 
[2OOII), which show a reasonable agreement with the measure- 
ments from the simulation with the exception of the high 
mass sample. Table [4] also shows the values of the Kaiser 
boost factor of equation p2p computed using h^s and the 
ones obtained from the measured power spectra (equation 
I14p : as for the dark matter case, they agree remarkably well. 
Since the dependency of equation p2|l on b is very weak, 
the difference between Sun computed with fceff and with b is 
much smaller than the error bars of Sat- The upper dotted 
lines in Figure [TT] correspond to the value of b^ff'S'iin- 



5 CONCLUSIONS 

The increasing quantity and quality of information from 
large galaxy redshift surveys demands models able to de- 
scribe the clustering of the galaxy distribution with high 
accuracy. The shape of the power spectrum, the tool most 
commonly used to analyze galaxy clustering, is distorted 
by non-linear evolution, redshift-space distortions and bias. 
These effects complicate the relation between the observa- 
tions and the predictions of linear perturbation theory, mak- 
ing the interpretation of these measurements in terms of 
constraints on cosmological parameters more difficult. 

In this paper we presented a model for the full 
shape of the large-scale power spectrum which is based 
on RPT and is analogous to t he approach fo l lowed by 
ICrocce fc Scoccimarrd (|2008l ) and ISanchez et al] (|2008l ) for 
the two-point correlation function. We compared our model 
against power spectra measured in a suite of 50 large vol- 
ume, moder ate resolution N-body simulations, call ed L- 
BASICC II (|Angulo et al.ll2008l : [SAnchez et al.ll2008l ). Our 
results indicate that the simple model presented here can 
provide an accurate description of the full shape of the 
power spectrum taking into account the effects of non-linear 
evolution, redshift-space distortions and halo bias for scales 
k < 0.15 /iMpc~^, making it a valuable tool for the analysis 
of real datasets. 

When applied to the dark matter distribution, our 
model gives a better description of the shape of the power 
spectrum at increasing redshift, where non- linear effects be- 
come less important. Even though this holds both in real 
and redshift-space, in the latter case the model performs 
worse, since it does not include explicitly the effect of scale 
dependent redshift-space distortions. The model also gives a 
correct description of the shape of the halo power spectrum 
for different mass ranges both for real and redshift-space 
information. In the latter case, due to the absence of the 
non-linear contribution from the fingers-of-god effect, the 
scale dependence of the redshift-space distortions is simpli- 
fied and the obtained constraints on the stretch parameter 
are similar to the ones of the real-space case. 

I Sanchez et al.l l|2008r ) performed a similar analysis for 
the shape of the large scale two-point correlation function. 
Their results showed that non-linear evolution, redshift- 
space distortions and bias are much more simpler to deal 
with for the case of the correlation function than for the 
power spectrum, where the signature is highly scale depen- 



dent. However, it is necessary to pursue complementary ap- 
proaches to constrain cosmological parameters allowing a 
comparison of the obtained results and a check for possi- 
ble systematics. The model for the shape of power spec- 
trum presented here has been proven to provide constraints 
on the stretch param eter a similar to the ones obtained by 
ISanchez et al.l (2008). Thus, with an accurate modelling, it is 
possible to extract cosmological information with the same 
precision from both statistics. 

Using a similar set of simulations to the ones used in 
this work. lAngulo et al.l (2008) tested a model in which the 
information from the broad band shape of the power spec- 
trum is discarded in order to extract a measurement of the 
BAO oscillations only. A comparison of their results with 
ours shows that the extra information in the full shape of the 
power spectrum helps to improve the obtained constraints 
over the BAO alone case (see also IShoii. Jeong fc Komatsrn 
l2009f ). In fact, despite the relatively small number of wave 
modes included in our analysis, the constrain ts that we ob- 
tain for a are slightly tighter than those of lAngulo et al.l 
(2008). The model equation (0 can be extended in a num- 
ber of ways, like for instance, with the use of the full RPT 
calculation of P2 mc and the inclusion of higher order terms 
in the description of Pmc or a specific treatment of the ef- 
fects of redshift-space distortions. This can help to extend 
the range of scales in which it is valid and improve the ac- 
curacy of the constraints on a and thus the dark energy 
equation of state. 
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APPENDIX A: POWER SPECTRUM: FFT AND 
CORRECTIONS 

The power spectrum is a statistical tool that allows to char- 
acterize gravitational clustering. In terms of the Fourier 
transform of the overdensity 5(x, i) of the continuous density 
field p{x,t), its definition is 



The Fourier transform of the MAS is given by 



(27r)^<5^(k - k)P(k,^) = {5{k, z)S*{L z)), 



(Al) 



where 5Q(k) is the three dimensional Dirac delta function. 
For a finite number A'^ of points tracing the underlying den- 
sity field in a limited volume V, the power spectrum becomes 
llJing||2005|): 



P(fc) = (|5(k)i^)- 



(A2) 



with n the mean number density of particles. The term 1/n 
is commonly called shot noise, is white, i.e. scale indepen- 
dent, and is given by the discretization of the density field. 
We Fourier transform the density field using a fast 
Fourier transform (FFT) algorithm that drastically in- 
creases the speed of the computation of the Fourier trans- 
form with respect to the standard transform. The drawback 
of the FFT is that it requires to convolve the density field 
5(x) with a regular grid. The conversion of 5(x) into a den- 
sity field 5'*(x) on a grid of cell size H is governed by the 
Mass Assignment Scheme (MAS) function VF(x). When the 
power spectrum is computed co nvolving the density with th e 
MAS, equation (|X2)) becomes (|Hocknev fc Eastwoodlll98ll ') 



|(5''(k)|^) = Y^ \W{]^ + 2fcNn)i^P(k + 2fcNn) 

n 

+ iy |w'(k+2fcNn)l^ 



(A3) 



where n is a three dimensional integer vector and ^n = ti / H 
is the Nyquist wavenumber and represent the larger value 
of k that is possible to achieve with the FFT. H^(k) is the 
Fourier transform of the MAS. 

Three MAS most commonly used are the Nearest 
Grid Point (NGP), Cloud in CeU (CIC) and Triangular 
Shape Cloud (TSC). In three dimensions, they are given 
by W(x) = nLi W{x,), where: 



WNGp(a;i) = 



Wcic{xi) 



1 |a^,|<0.5, 
\x,\ >0.5, 



<1, 
> 1, 



WTSc(a;i) = 



0.75 -a;f \xi\<Q.^, 
iL^^hMl 0.5 s; \x,\ < 1.5, 
,0 bd^l.S. 



(A4a) 



(A4b) 



(A4c 



W"(k) 



n 



i(7rfci/2feN) 



7^fci/2fc^ 



(A4d) 



with p the order of the MAS (p = 1 for NCC, p = 2 for CIC 
andp = 3for TSC). 

The computation of the power spectrum with a FFT 
algorithm is divided roughly in the following steps : 

(i) creation of a density field assigning the particles to a 
grid using a chosen MAS; 

(ii) Fourier transform of the density field - we have per- 
formed this step using the free software FFTW; 

(iii) spherical average of the density field in Fourier space 
to obtain the left hand side of equation (|A3[) : 

(iv) solve the equation (|A3|) to obtain P{k). 

In Ijingj (120051 ). the last point is solved with an itera- 
tive procedure based on the assumption that the PS can be 
approximated by a power law at fc > fcN. This assumption 
is a viable approximation when computing the power spec- 
trum for dark matter, but it is not applicable when dealing 
with the halo-halo power spectrum (see Figure [2] and ii2.2|) . 
It is also possible to approxima te the solution dividing the 
computed {5'*(k)) by VF(k) ( e.g.lAngulo et al.l 12008)) or by 
(E„|l^(k + 2fcNn)|2)i/2 (e.g lJeong fc Komatsullioog )). So 
the last two steps of the computation become: 

(iii) correction of the amplitude of each Fourier mode for 
the effect of the MASS; 

(iv) spherical average of the density field in Fourier space 
to obtain the left hand side of equation (|X3]). 



The recovered power spectra for the two corrections de- 
scribed above are 



{\5\\.)\^) 1 
^'^^>- W2(k) n-^^^> 






and 



P2(fc) = 



■|VK(k') 


"P(k') 

(k) 


1 

+ - 
n 


M/(k')in 

Vl/2(k) _ 


{|5'(k) 


') 


1 


r>n.\ 



(A5a 



X:„lVK(k + 2fcNn)|2 n 



1 + 



E|„|^olW^(k')Pnk') 



W^{k)P{k) 
EH^o!l^(k')" 



(A5b) 



W2(k) 



+ 



where k' = k -f 2fcNn. Equation ()A5b|) makes use of the 
fact that En!''^(k')P is finite and that for jkj < /cn, 
E|„|^o |W^(k')l < |W^(k)i and P(k') ^ P(k). 

We have tested the impact of the MAS and correc- 
tions on the recovered power spectrum. To do this we have 
compared the power spectrum computed with the standard 
Fourier transform (FT) and with the FFTW algorithm. For 
the FFT we used a grid of 200^ cells, filled using one of 
the three possible MAS. We have also switched on and off 
the corrections in equations (|A5|I . In order to have the same 
number of modes in the two cases, the values of k used for 
the FT calculation are the same as for the FFT. The power 
spectrum for all the cases has been computed in 48 logarith- 
mic bins between k = 0.01 ZiMpc"^ and fcN = 0.47 /iMpc^^. 
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Figure Al. Ratio between the power spectrum computed using tiie FFT and tiie one computed with the standard FT versus the 
wavenumber k. For clearness the error bars are not shown (for fc > 0.03/iMpc~^ their amphtude is comparable with the thickness of 
the lines). The FFT has been computed on a grid of 200^ cells, that has been filled using the NGP (left panel), CIC (central panel) and 
TSC (right panel). The solid black line is the ratio between the uncorrected PppT(fc) and the PpT(fc), while the dashed red line and the 
dash-dotted blue line indicate the Ppft(A:) corrected as in equations l)A5a|l and l|A5b| l respectively. Since J^^ |iy(k + 2A:Nn)p = 1, P2{k) 
is equivalent to the uncorrected power spectrum in the NGP case. The dotted vertical lino outlines the 67%feN. 

Figure lAll shows the results of this test as the ratio be- 
tween the FFT and the FT power spectra as function of the 
wave-number k. The panels a, b and c correspond, respec- 
tively, to the power spectrum computed with NGP, CIC 
and TSC as MAS. In each panel the solid line is for the 
uncorrected power spectrum Punc(fc) = (|'5'*(fc)P) — 1/n, 
the dashed line is for the Pi{k) while the dash-dotted one 
is for P2{k). Since J2n 1^0^ + 2fcNn)p = 1 in the NGP 
case, P2{k) = Punc(fc). The correction of the {\5^{k)\^) in- 
creases drastically the agreement with the real P{k), and 
the achieved accuracy is larger the higher the order of the 
MAS is; equation (IA5b[l seems to work slightly better than 
equation (|A5a|l . In the light of these results, we have decided 
to use the TSC as MAS correcting the power spectrum as in 
equation (|A5b|l for our analysis. Choosing a very conserva- 
tive limit, we assume that the FFT gives a correct answer for 
the power spectrum if PFFT{k)/PFT{k) < 0.5%: this holds 
for k ~ 0.314 /iMpc^^ — 67%kN, which is indicated by the 
vertical dotted line in Figure IXT] 

This paper has been typeset from a T^jX/ I^TJ^X file prepared 
by the author. 
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